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^ ; We perform a NLO QCD analysis of the combined SLAC-BCDMS-NMC-E665-H1-ZEUS data 

on inclusive deep inelastic cross section. Particular attention was paid to the extraction of 

strong coupling constant Os and high twist (HT) contribution to the structure functions F2 and 

^ , -Fl . It was shown that at small and moderate x there is a visible dependence of the extracted 

^^ ■ values of HT contribution to F2 on the QCD renormalization scale, which indicates that in 

Q/ ' this region extracted HT can absorb NNLO QCD corrections. At larger x the dependence of 

HT on the renormalization scale is negligible and the influence of NNLO correction on their 

values should be less significant. The value of as{Mz) = 0.1159 ± 0.0031 (total) is obtained, 

where the error includes statistical, systematical and theoretical uncertainties. 
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Q.' The problem of strong coupling constant determination from charged leptons deep inelastic 

scattering (DIS) data was widely discussed recent years. These data are very precise (at the 

Q I level of 0(1%)) and the theoretical uncertainties of the analysis are relatively small, which 

i-G ' allows for to determine the value of as{Mz) with the precision of 0(0.001). At the same time 

there is statistically significant discrepancy between the value of a^ obtained in the analysis H of 
k><( \ combined SLAC-BCDMS proton-deuterium dataH'a and the results of experiments performed 

\^ ' at LEPq. This discrepancy can be considered as an indication on new physics beyond Standard 

ModelEl. Meanwhile, as it became clear in the earliest QCD analysis of DIS data the value of Ug 
is strongly correlated with the value of possible high twist (HT) contribution to the structure 
function F2 H. This correlation makes the separation of log-like and power-like contributions 
to the scaling violation unstable with respect to various assumptions made in the analysis. In 
particular, as it was shown recentlyQ, the results of nonsinglet SLAC-BCDMS proton-deuterium 
data analysis are sensitive to the procedures used to handle systematic errors on the data. 
The central value of as{Mz) = 0.1180 ±0.0017 (stat.+syst.), as obtained in the analysis of RefQ 
with the complete account ofpoint-to-point correlations due to systematic errors, is significantly 
larger than the results of RefH and is compatible with the LEP measurements and world average. 
In the extended version of this analysis with addition of the NMC proton-deuterium dataa and 
account of HT contribution to the structure function Fl the value of as{Mz) = 0.1170 ± 0.0021 
(stat.+syst.) was obtained H. In this talk we describe the effect of further extension of the 
analysed data set on the value of as{Mz)- 

The analysis was performed in NLO QCD approximation with fixed number of evolved 
fermion distributions (A'^^ = 3); momentum sum rule (MSR) and fermion sum rule for valence 
quarks were used to decrease the number of fitted parameters; HT contributions to the structure 
functions F2 and F^ were parametrized in additive form; target mass correctionll2l was accounted 
for in the fitted cross section formula; account of systematic errors was performed through the 



covariance matrix approach. More detailed description of the ansatz can be found in the earher 
papers Qtl'Ell. The milestone results are given in Table 1. From the top to bottom rows the 
precision and reliability of the a^ improves due to more data included in the analysis and more 
corrections applied. The cut Q^ > 2.5 GeV^ was imposed because without this cut the total 
error is dominated by the uncertainty in QCD renormalization scale (see fourth row of Table 1) 
and it is meaningless to suppress the experimental error with a risk to encounter an unexpected 
small Q^ effects. Our final value of as{Mz) is given in the last row of Table 1 and its various 
theoretical uncertainties - in Table 2. With the account of the shifts due to this uncertainties 
one can obtain 

as{Mz) = 0.1159 ±0.0031(total), 

that is compatible with world average. In the analysis of Ref§ it was obtained, that in the 
QCD fit to the world DIS data, performed without low Q^ cut, MSR is violated. In order to 
check this conclusion a test fit without imposing MSR boundary condition was made. Obtained 
value of total momentum carried by partons is < x >= 0.982 ± 0.028 that is compatible with 1 
within errors and is in disagreement with the results of RefO (< x >~ 1.08 ±0.02 for the cut 
Q2 > 3 GeV2). 



Table 1: The values of as{Mz) obtained in the analysis of different data sets. NDP is total number of data points, 
experimental error (exp) includes statistical and systematical uncertainties, (RS) is the shift due to the change of 

QCD renormalization scale from Q^ to 4Q^. 



Analysed data set 



SLAC-BCDMS-NMC (0.3 < x < 0.75) 



SLAC-BCDMS-NMC (0.3 <x< 0.75, 
Fermi motion correction 113 (FMC) on) 



SLAC-BCDMS-NMC (x > 0.3, FMC on) 



SLAC-BCDMS-NMC (no x-cut, FMC on) 



SLAC-BCDMS-NMC 
(no x-cut, Q2 > 2.5 GeYi, FI 



sVl, FMC 



.Con) 



SLAC-BCDMS-NMC-E66i^-H£l-ZEU 
no x-cut, Q2 > 2.5 GeV^, FMC on 



# 



NDP 



1243 



1243 



1348 



2541 



2083 



2512 



as{Mz 



0.1170 ±0.0021(exp) 



0.1190 ±0.0020(exp) 



0.1197 ±0.0019(exp) 



0.1190 ±0.0012(exp)± 
+0.0028(RS) 



0.1170 ± 0.0019(exp) 



0.1166 ±0.0016 (exp) 



Table 2: The effect of various theoretical uncertainties on the value of as{Mz) 



Source of uncertainty 


AasiMz) 


QCD renormalization scale variation from l/AQ'^ to 4Q"^ 


-0.0022, ±0.0028 


The change of matching scale of heavy 
quark threshold from mq to 6.5mQ(c.f. RefUB) 


-0.0020 


The change of c-quark mass on ±0.25 GeV 


±0.0002 


The change of strange sea suppression factor on ±0.1 


±0.0001 



The fitted values of HT contributions to the proton structure function F2 and nucleon 
structure function Fl for different choices of the QCD renormalization scale // are given in Fig. 
1. (It was adopted in the analysis that HT contributions to the proton and neutron structure 
functions Fl are equal, since data cannot discriminate between them). The visible dependence 
of H2{x) on /x at small x can be considered as an indication that in this region HT contribution 
to F2 absorbs NNLO QCD corrections. The indication on interplay between HT contribution 
to the structure function F^ and NNLO corrections was obtained in the NLO QCD analysis 
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Figure 1: The fitted values of Hf (x) and H^ (x) for different choices of the QCD renormalization scale fi. Full 
circles correspond to /x = Q^, empty circles - to /i = 4Q^, squares - to jU = Q^/4 

of neutrino data t3; this interplay was also directly demonstrated in the earlier NNLO QCD 
analysis of RefEj. At the same time for H2{x) at largest x and for Hl{x) this dependence is 
not so significant. It is worth to note that in the fit with simultaneous extraction of a^ and HT 
contribution the dependence of the Og value on /i is weaker, than in the fit with HT fixed since 
HT are readjusted with the change of fi; see in this connection Fig. 3 of RefB3. In particular, 
dueJio this effect our value of the RS uncertainty on as{Mz) is less than it was obtained in 
RefE-3. At the same time due to the large correlation between ag and HT contribution to F2 the 
as error increases as compared to the fit with HT fixed and hence one can say that some part 
of the RS error on ct^ obtained in the fit with HT released is included in the total experimental 



error. 



I am indebted to G.Korchemsky and L.Mankiewicz for useful discussions. 
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